The Puerto Rican Infant Metagenomic and Epidemiologic Study of Respiratory Outcomes (PRIMERO): Design and Baseline Characteristics for a Birth Cohort Study of Early-life Viral Respiratory Illnesses and Airway Dysfunction in Puerto Rican Children.

Epidemiologic studies demonstrate an association between early-life respiratory illnesses (RIs) and the development of childhood asthma. However, it remains uncertain whether these children are predisposed to both conditions or if early-life RIs induce alterations in airway function, immune responses, or other human biology that contribute to the development of asthma. Puerto Rican children experience a disproportionate burden of early-life RIs and asthma, making them an important population for investigating this complex interplay. PRIMERO, the Puerto Rican Infant Metagenomics and Epidemiologic Study of Respiratory Outcomes , recruited pregnant women and their newborns to investigate how the airways develop in early life among infants exposed to different viral RIs, and will thus provide a critical understanding of childhood asthma development. As the first asthma birth cohort in Puerto Rico, PRIMERO will prospectively follow 2,100 term healthy infants. Collected samples include post-term maternal peripheral blood, infant cord blood, the child's peripheral blood at the year two visit, and the child's nasal airway epithelium, collected using minimally invasive nasal swabs, at birth, during RIs over the first two years of life, and at annual healthy visits until age five. Herein, we describe the study's design, population, recruitment strategy, study visits and procedures, and primary outcomes.

A common polymorphism in the Intelectin-1 gene influences mucus plugging in severe asthma.

By incompletely understood mechanisms, type 2 (T2) inflammation present in the airways of severe asthmatics drives the formation of pathologic mucus which leads to airway mucus plugging. Here we investigate the molecular role and clinical significance of intelectin-1 (ITLN-1) in the development of pathologic airway mucus in asthma. Through analyses of human airway epithelial cells we find that ITLN1 gene expression is highly induced by interleukin-13 (IL-13) in a subset of metaplastic MUC5AC+ mucus secretory cells, and that ITLN-1 protein is a secreted component of IL-13-induced mucus. Additionally, we find ITLN-1 protein binds the C-terminus of the MUC5AC mucin and that its deletion in airway epithelial cells partially reverses IL-13-induced mucostasis. Through analysis of nasal airway epithelial brushings, we find that ITLN1 is highly expressed in T2-high asthmatics, when compared to T2-low children. Furthermore, we demonstrate that both ITLN-1 gene expression and protein levels are significantly reduced by a common genetic variant that is associated with protection from the formation of mucus plugs in T2-high asthma. This work identifies an important biomarker and targetable pathways for the treatment of mucus obstruction in asthma.

Atopic dermatitis complicated by recurrent eczema herpeticum is characterized by multiple, concurrent epidermal inflammatory endotypes.

BACKGROUND: A subgroup of atopic dermatitis (AD) patients suffer from recurrent, disseminated herpes simplex virus (HSV) skin infections, termed eczema herpeticum (EH), which can be life-threatening and contribute to AD morbidity. The pathobiology underlying ADEH is unknown. OBJECTIVE: To determine transcriptional mechanisms of skin and immune system pathobiology that underlie ADEH disease. METHODS: We performed whole transcriptome RNA-sequencing of non-lesional skin samples (epidermis, dermis) of AD patients with (ADEH + , n=15) and without (ADEH - , n=13) recurrent EH history, and healthy controls (HC, n=15). We also performed RNA-sequencing on plasmacytoid dendritic cells (pDCs) collected from these participants and infected in vitro with HSV-1. Differential expression, gene set enrichment, and endotyping analyses were performed. RESULTS: ADEH + disease was characterized by dysregulation in skin gene expression, which was limited in dermis (differentially expressed genes [DEGs]=14) and widespread in epidermis (DEGs=129). ADEH + -upregulated epidermal DEGs were enriched in type 2 cytokine (T2) ( IL4R, CCL22, CRLF2, IL7R ), interferon ( CXCL10, ICAM1, IFI44 , and IRF7) , and IL-36γ ( IL36G ) inflammatory pathway genes. At a person-level, all ADEH + participants exhibited T2 and interferon endotypes and 87% were IL36G-high. In contrast, these endotypes were more variably expressed among ADEH - participants. ADEH + patient skin also exhibited dysregulation in epidermal differentiation complex (EDC) genes within the LCE, S100 , and SPRR families, which are involved in skin barrier function, inflammation, and antimicrobial activities. pDC transcriptional responses to HSV-1 infection were not altered by ADEH status. CONCLUSIONS: ADEH + pathobiology is characterized by a unique, multi-faceted epidermal inflammation that accompanies dysregulation in the expression of EDC genes. Key Messages: AD patients with a history of recurrent EH exhibit molecular skin pathobiology that is similar in form, but more severe in degree, than in AD patients without this complication. Non-lesional skin of ADEH + patients concurrently exhibits excessive type 2 cytokine, interferon, and IL-36γ-driven epidermal inflammation. Expression of these inflammatory skin endotypes among ADEH + patients is associated with dysregulation in expression of epidermal differentiation complex genes involved in barrier function, inflammation, and antimicrobial activity. Capsule Summary: AD patients with a history of recurrent disseminated HSV-1 skin infections form a unique molecular skin endotype group that concurrently exhibits type 2 cytokine, interferon, and IL-36γ-driven skin inflammation, accompanied by dysregulation in expression of epidermal differentiation complex genes involved in barrier function, inflammation, and antimicrobial activity.

Nasal airway transcriptome-wide association study of asthma reveals genetically driven mucus pathobiology.

To identify genetic determinants of airway dysfunction, we performed a transcriptome-wide association study for asthma by combining RNA-seq data from the nasal airway epithelium of 681 children, with UK Biobank genetic association data. Our airway analysis identified 95 asthma genes, 58 of which were not identified by transcriptome-wide association analyses using other asthma-relevant tissues. Among these genes were MUC5AC, an airway mucin, and FOXA3, a transcriptional driver of mucus metaplasia. Muco-ciliary epithelial cultures from genotyped donors revealed that the MUC5AC risk variant increases MUC5AC protein secretion and mucus secretory cell frequency. Airway transcriptome-wide association analyses for mucus production and chronic cough also identified MUC5AC. These cis-expression variants were associated with trans effects on expression; the MUC5AC variant was associated with upregulation of non-inflammatory mucus secretory network genes, while the FOXA3 variant was associated with upregulation of type-2 inflammation-induced mucus-metaplasia pathway genes. Our results reveal genetic mechanisms of airway mucus pathobiology.

Single cell analysis of host response to helminth infection reveals the clonal breadth, heterogeneity, and tissue-specific programming of the responding CD4+ T cell repertoire.

The CD4+ T cell response is critical to host protection against helminth infection. How this response varies across different hosts and tissues remains an important gap in our understanding. Using IL-4-reporter mice to identify responding CD4+ T cells to Nippostrongylus brasiliensis infection, T cell receptor sequencing paired with novel clustering algorithms revealed a broadly reactive and clonally diverse CD4+ T cell response. While the most prevalent clones and clonotypes exhibited some tissue selectivity, most were observed to reside in both the lung and lung-draining lymph nodes. Antigen-reactivity of the broader repertoires was predicted to be shared across both tissues and individual mice. Transcriptome, trajectory, and chromatin accessibility analysis of lung and lymph-node repertoires revealed three unique but related populations of responding IL-4+ CD4+ T cells consistent with T follicular helper, T helper 2, and a transitional population sharing similarity with both populations. The shared antigen reactivity of lymph node and lung repertoires combined with the adoption of tissue-specific gene programs allows for the pairing of cellular and humoral responses critical to the orchestration of anti-helminth immunity.

Type 2 and interferon inflammation regulate SARS-CoV-2 entry factor expression in the airway epithelium.

Coronavirus disease 2019 (COVID-19) is caused by SARS-CoV-2, an emerging virus that utilizes host proteins ACE2 and TMPRSS2 as entry factors. Understanding the factors affecting the pattern and levels of expression of these genes is important for deeper understanding of SARS-CoV-2 tropism and pathogenesis. Here we explore the role of genetics and co-expression networks in regulating these genes in the airway, through the analysis of nasal airway transcriptome data from 695 children. We identify expression quantitative trait loci for both ACE2 and TMPRSS2, that vary in frequency across world populations. We find TMPRSS2 is part of a mucus secretory network, highly upregulated by type 2 (T2) inflammation through the action of interleukin-13, and that the interferon response to respiratory viruses highly upregulates ACE2 expression. IL-13 and virus infection mediated effects on ACE2 expression were also observed at the protein level in the airway epithelium. Finally, we define airway responses to common coronavirus infections in children, finding that these infections generate host responses similar to other viral species, including upregulation of IL6 and ACE2. Our results reveal possible mechanisms influencing SARS-CoV-2 infectivity and COVID-19 clinical outcomes.

Single-Cell and Population Transcriptomics Reveal Pan-epithelial Remodeling in Type 2-High Asthma.

The type 2 cytokine-high asthma endotype (T2H) is characterized by IL-13-driven mucus obstruction of the airways. To further investigate this incompletely understood pathobiology, we characterize IL-13 effects on human airway epithelial cell cultures using single-cell RNA sequencing, finding that IL-13 generates a distinctive transcriptional state for each cell type. Specifically, we discover a mucus secretory program induced by IL-13 in all cell types which converts both mucus and defense secretory cells into a metaplastic state with emergent mucin production and secretion, while leading to ER stress and cell death in ciliated cells. The IL-13-remodeled epithelium secretes a pathologic, mucin-imbalanced, and innate immunity-depleted proteome that arrests mucociliary motion. Signatures of IL-13-induced cellular remodeling are mirrored by transcriptional signatures characteristic of the nasal airway epithelium within T2H versus T2-low asthmatic children. Our results reveal the epithelium-wide scope of T2H asthma and present candidate therapeutic targets for restoring normal epithelial function.

Type 2 and interferon inflammation strongly regulate SARS-CoV-2 related gene expression in the airway epithelium.

Coronavirus disease 2019 (COVID-19) outcomes vary from asymptomatic infection to death. This disparity may reflect different airway levels of the SARS-CoV-2 receptor, ACE2, and the spike protein activator, TMPRSS2. Here we explore the role of genetics and co-expression networks in regulating these genes in the airway, through the analysis of nasal airway transcriptome data from 695 children. We identify expression quantitative trait loci (eQTL) for both ACE2 and TMPRSS2, that vary in frequency across world populations. Importantly, we find TMPRSS2 is part of a mucus secretory network, highly upregulated by T2 inflammation through the action of interleukin-13, and that interferon response to respiratory viruses highly upregulates ACE2 expression. Finally, we define airway responses to coronavirus infections in children, finding that these infections upregulate IL6 while also stimulating a more pronounced cytotoxic immune response relative to other respiratory viruses. Our results reveal mechanisms likely influencing SARS-CoV-2 infectivity and COVID-19 clinical outcomes.

Dissecting the cellular specificity of smoking effects and reconstructing lineages in the human airway epithelium.

Cigarette smoke first interacts with the lung through the cellularly diverse airway epithelium and goes on to drive development of most chronic lung diseases. Here, through single cell RNA-sequencing analysis of the tracheal epithelium from smokers and non-smokers, we generate a comprehensive atlas of epithelial cell types and states, connect these into lineages, and define cell-specific responses to smoking. Our analysis infers multi-state lineages that develop into surface mucus secretory and ciliated cells and then contrasts these to the unique specification of submucosal gland (SMG) cells. Accompanying knockout studies reveal that tuft-like cells are the likely progenitor of both pulmonary neuroendocrine cells and CFTR-rich ionocytes. Our smoking analysis finds that all cell types, including protected stem and SMG populations, are affected by smoking through both pan-epithelial smoking response networks and hundreds of cell-specific response genes, redefining the penetrance and cellular specificity of smoking effects on the human airway epithelium.

COVID-19-related Genes in Sputum Cells in Asthma. Relationship to Demographic Features and Corticosteroids.

Rationale: Coronavirus disease (COVID-19) is caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). ACE2 (angiotensin-converting enzyme 2), and TMPRSS2 (transmembrane protease serine 2) mediate viral infection of host cells. We reasoned that differences in ACE2 or TMPRSS2 gene expression in sputum cells among patients with asthma may identify subgroups at risk for COVID-19 morbidity.Objectives: To determine the relationship between demographic features and sputum ACE2 and TMPRSS2 gene expression in asthma.Methods: We analyzed gene expression for ACE2 and TMPRSS2, and for ICAM-1 (intercellular adhesion molecule 1) (rhinovirus receptor as a comparator) in sputum cells from 330 participants in SARP-3 (Severe Asthma Research Program-3) and 79 healthy control subjects.Measurements and Main Results: Gene expression of ACE2 was lower than TMPRSS2, and expression levels of both genes were similar in asthma and health. Among patients with asthma, male sex, African American race, and history of diabetes mellitus were associated with higher expression of ACE2 and TMPRSS2. Use of inhaled corticosteroids (ICS) was associated with lower expression of ACE2 and TMPRSS2, but treatment with triamcinolone acetonide did not decrease expression of either gene. These findings differed from those for ICAM-1, where gene expression was increased in asthma and less consistent differences were observed related to sex, race, and use of ICS.Conclusions: Higher expression of ACE2 and TMPRSS2 in males, African Americans, and patients with diabetes mellitus provides rationale for monitoring these asthma subgroups for poor COVID-19 outcomes. The lower expression of ACE2 and TMPRSS2 with ICS use warrants prospective study of ICS use as a predictor of decreased susceptibility to SARS-CoV-2 infection and decreased COVID-19 morbidity.

Genome-Wide Analysis Reveals Mucociliary Remodeling of the Nasal Airway Epithelium Induced by Urban PM2.5.

Air pollution particulate matter <2.5 µm (PM2.5) exposure is associated with poor respiratory outcomes. Mechanisms underlying PM2.5-induced lung pathobiology are poorly understood but likely involve cellular and molecular changes to the airway epithelium. We extracted and chemically characterized the organic and water-soluble components of air pollution PM2.5 samples, then determined the whole transcriptome response of human nasal mucociliary airway epithelial cultures to a dose series of PM2.5 extracts. We found that PM2.5 organic extract (OE), but not water-soluble extract, elicited a potent, dose-dependent transcriptomic response from the mucociliary epithelium. Exposure to a moderate OE dose modified the expression of 424 genes, including activation of aryl hydrocarbon receptor signaling and an IL-1 inflammatory program. We generated an OE-response gene network defined by eight functional enrichment groups, which exhibited high connectivity through CYP1A1, IL1A, and IL1B. This OE exposure also robustly activated a mucus secretory expression program (>100 genes), which included transcriptional drivers of mucus metaplasia (SPDEF and FOXA3). Exposure to a higher OE dose modified the expression of 1,240 genes and further exacerbated expression responses observed at the moderate dose, including the mucus secretory program. Moreover, the higher OE dose significantly increased the MUC5AC/MUC5B gel-forming mucin expression ratio and strongly downregulated ciliated cell expression programs, including key ciliating cell transcription factors (e.g., FOXJ1 and MCIDAS). Chronic OE stimulation induced mucus metaplasia-like remodeling characterized by increases in MUC5AC+ secretory cells and MUC5AC mucus secretions. This epithelial remodeling may underlie poor respiratory outcomes associated with high PM2.5 exposure.

Functional genomics of CDHR3 confirms its role in HRV-C infection and childhood asthma exacerbations.

BACKGROUND: Research in transformed immortalized cell lines indicates the cadherin-related family member 3 (CDHR3) protein serves as a receptor for human rhinovirus (HRV)-C. Similar experiments indicate that the CDHR3 coding variant rs6967330 increases CDHR3 protein surface expression. OBJECTIVE: We sought to determine whether CDHR3 is necessary for HRV-C infection of primary airway epithelial cells (AECs) and to identify molecular mechanisms by which CDHR3 variants confer risk for asthma exacerbations. METHODS: CDHR3 function and influence on HRV-C infection were investigated by using single-cell transcriptomics, CRISPR-Cas9 gene knockout, and genotype-specific donor experiments performed in primary AECs. Nasal airway epithelium cis-expression quantitative trait locus (eQTL) analysis of CDHR3 was performed, followed by association testing for asthma hospitalization in minority children. RESULTS: CDHR3 lung expression is exclusive to ciliated AECs and associated with basal bodies during and after motile ciliogenesis. Knockout of CDHR3 in human AECs did not prevent ciliated cell differentiation but was associated with a decrease in transepithelial resistance and an 80% decrease in HRV-C infection of the mucociliary epithelium. AECs from subjects homozygous for the risk-associated rs6967330 single nucleotide polymorphism (SNP) exhibited greater HRV-C infection compared with cells homozygous for the nonrisk allele. AEC cis-eQTL analysis indicated that rs6967330 and other SNPs are eQTLs for CDHR3. Only the eQTL block containing the rs6967330 SNP showed a significant association with childhood asthma hospitalization. CONCLUSIONS: Genetic deletion and genotype-specific studies in primary AECs indicate CDHR3 is critical to HRV-C infection of ciliated cells. The rs6967330 SNP confers risk of severe childhood asthma exacerbations, likely through increasing HRV-C infection levels and protein surface localization.

The nonlesional skin surface distinguishes atopic dermatitis with food allergy as a unique endotype.

Skin barrier dysfunction has been reported in both atopic dermatitis (AD) and food allergy (FA). However, only one-third of patients with AD have FA. The purpose of this study was to use a minimally invasive skin tape strip sampling method and a multiomics approach to determine whether children with AD and FA (AD FA+) have stratum corneum (SC) abnormalities that distinguish them from AD without FA (AD FA-) and nonatopic (NA) controls. Transepidermal water loss was found to be increased in AD FA+. Filaggrin and the proportion of ω-hydroxy fatty acid sphingosine ceramide content in nonlesional skin of children with AD FA+ were substantially lower than in AD FA- and NA skin. These abnormalities correlated with morphologic changes in epidermal lamellar bilayer architecture responsible for barrier homeostasis. Shotgun metagenomic studies revealed that the nonlesional skin of AD FA+ had increased abundance of Staphylococcus aureus compared to NA. Increased expression of keratins 5, 14, and 16 indicative of hyperproliferative keratinocytes was observed in the SC of AD FA+. The skin transcriptome of AD FA+ had increased gene expression for dendritic cells and type 2 immune pathways. A network analysis revealed keratins 5, 14, and 16 were positively correlated with AD FA+, whereas filaggrin breakdown products were negatively correlated with AD FA+. These data suggest that the most superficial compartment of nonlesional skin in AD FA+ has unique properties associated with an immature skin barrier and type 2 immune activation.

CD11c+ Cells Are Gatekeepers for Lymphocyte Trafficking to Infiltrated Islets During Type 1 Diabetes.

Type 1 diabetes (T1D) is a T cell mediated autoimmune disease that affects more than 19 million people with incidence increasing rapidly worldwide. For T cells to effectively drive T1D, they must first traffic to the islets and extravasate through the islet vasculature. Understanding the cues that lead to T cell entry into inflamed islets is important because diagnosed T1D patients already have established immune infiltration of their islets. Here we show that CD11c+ cells are a key mediator of T cell trafficking to infiltrated islets in non-obese diabetic (NOD) mice. Using intravital 2-photon islet imaging we show that T cell extravasation into the islets is an extended process, with T cells arresting in the islet vasculature in close proximity to perivascular CD11c+ cells. Antigen is not required for T cell trafficking to infiltrated islets, but T cell chemokine receptor signaling is necessary. Using RNAseq, we show that islet CD11c+ cells express over 20 different chemokines that bind chemokine receptors expressed on islet T cells. One highly expressed chemokine-receptor pair is CXCL16-CXCR6. However, NOD. CXCR6-/- mice progressed normally to T1D and CXCR6 deficient T cells trafficked normally to the islets. Even with CXCR3 and CXCR6 dual deficiency, T cells trafficked to infiltrated islets. These data reinforce that chemokine receptor signaling is highly redundant for T cell trafficking to inflamed islets. Importantly, depletion of CD11c+ cells strongly inhibited T cell trafficking to infiltrated islets of NOD mice. We suggest that targeted depletion of CD11c+ cells associated with the islet vasculature may yield a therapeutic target to inhibit T cell trafficking to inflamed islets to prevent progression of T1D.

Expression and function of the ectopic olfactory receptor OR10G7 in patients with atopic dermatitis.

BACKGROUND: Ectopic olfactory receptors (ORs) are found in the skin, but their expression and biological function in normal skin and skin form patients with atopic dermatitis (AD) are unknown. OBJECTIVES: We sought to characterize the expression of ORs in the skin and assess OR-mediated biological responses of primary human keratinocytes in the presence of odorant ligands. METHODS: OR expression was examined by using whole-transcriptome sequencing of skin tape strips collected from patients with AD and healthy control (HC) subjects. OR10G7 and filaggrin 1 (FLG-1) expression was analyzed by using RT-PCR and immunostaining in skin biopsy specimens and primary human keratinocytes from patients with AD and HC subjects. ATP and cyclic AMP production by control and OR10G7 small interfering RNA-transfected keratinocytes in response to odorant stimulation with acetophenone and eugenol was assessed. RESULTS: A total of 381 OR gene transcripts were detected in the skin samples, with the greatest OR expression detected in the skin tape strips corresponding to the upper granular layer of the skin. OR10G7 expression was significantly increased in skin biopsy specimens from patients with AD compared with those from HC subjects (P = .01) and inversely correlated with FLG-1 expression (P = .009). OR10G7 expression was greatest in undifferentiated keratinocytes from patients with AD and was downregulated with progressive differentiation. Primary human keratinocytes produced ATP, an essential neurotransmitter in sensory pathways, in response to acetophenone and eugenol, odorants previously identified as potential ligands for this receptor. This response was abolished in OR10G7 small interfering RNA-transfected keratinocytes. CONCLUSIONS: OR10G7 is expressed at significantly greater levels in undifferentiated keratinocytes from patients with AD compared with HC subjects. OR10G7 is likely involved in transmission of skin-induced chemosensory responses to odorant stimulation, which might modulate differential nociceptive responses in AD skin.

A Transcriptomic Method to Determine Airway Immune Dysfunction in T2-High and T2-Low Asthma.

BACKGROUND: Type 2 (T2) inflammation drives airway dysfunction in many patients with asthma; yet, we lack a comprehensive understanding of the airway immune cell types and networks that sustain this inflammation. Moreover, defects in the airway immune system in patients with asthma without T2 inflammation are not established. OBJECTIVES: To determine the gene networks that sustain T2 airway inflammation in T2-high asthma and to explore the gene networks that characterize T2-low asthma. METHODS: Network analysis of sputum cell transcriptome expression data from 84 subjects with asthma and 27 healthy control subjects was used to identify immune cell type-enriched networks that underlie asthma subgroups. RESULTS: Sputum T2 gene expression was characterized by an immune cell network derived from multiple innate immune cells, including eosinophils, mast cells/basophils, and inflammatory dendritic cells. Clustering of subjects within this network stratified subjects into T2-high and T2-low groups, but it also revealed a subgroup of T2-high subjects with uniformly higher expression of the T2 network. These "T2-ultrahigh subjects" were characterized clinically by older age and more severe airflow obstruction and pathologically by a second T2 network derived from T2-skewed, CD11b+/CD103-/IRF4+ classical dendritic cells. Subjects with T2-low asthma were differentiated from healthy control subjects by lower expression of a cytotoxic CD8+ T-cell network, which was negatively correlated with body mass index and plasma IL-6 concentrations. CONCLUSIONS: Persistent airway T2 inflammation is a complex construct of innate and adaptive immunity gene expression networks that are variable across individuals with asthma and persist despite steroid treatment. Individuals with T2-low asthma exhibit an airway deficiency in cytotoxic T cells associated with obesity-driven inflammation.

Utilization of Air-Liquid Interface Cultures as an In Vitro Model to Assess Primary Airway Epithelial Cell Responses to the Type 2 Cytokine Interleukin-13.

The airway epithelium lines the respiratory tract and provides the primary protective barrier against inhalational insults including toxic environmental substances and microorganisms. The airway epithelium also plays a critical role in regulating airway immune responses. The airway epithelial response to the type 2 cytokine, interleukin-13 (IL-13), is critical to airway inflammation, mucus production, and airway hyperresponsiveness present in asthma. Relevant primary cell models of the human airway epithelium are needed to investigate the biology of IL-13-mediated airway epithelial effects. Here, we describe the generation of a differentiated mucociliary human airway epithelium using an in vitro air-liquid interface (ALI) culture model system. We also describe methods to stimulate this culture model with IL-13 and harvest cells and biomolecules to interrogate cellular and molecular aspects of the airway epithelial IL-13 response.

Correction to: Dual RNA-seq reveals viral infections in asthmatic children without respiratory illness which are associated with changes in the airway transcriptome.

In our recent article [1], it has come to our attention that the sample labels are not consistent between Table 1, the data labels deposited in the Sequence Read Archive, and Additional file 1: Table S2. We are therefore providing an updated Additional file 1: Table S2 so identical samples now have the same label.

ROP: dumpster diving in RNA-sequencing to find the source of 1 trillion reads across diverse adult human tissues.

High-throughput RNA-sequencing (RNA-seq) technologies provide an unprecedented opportunity to explore the individual transcriptome. Unmapped reads are a large and often overlooked output of standard RNA-seq analyses. Here, we present Read Origin Protocol (ROP), a tool for discovering the source of all reads originating from complex RNA molecules. We apply ROP to samples across 2630 individuals from 54 diverse human tissues. Our approach can account for 99.9% of 1 trillion reads of various read length. Additionally, we use ROP to investigate the functional mechanisms underlying connections between the immune system, microbiome, and disease. ROP is freely available at https://github.com/smangul1/rop/wiki .